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^ ■ I. Introduction 

Quantum dots (i.e. nanostructures with charge carriers confinement in all spatial directions) 
have an atom-like energy spectrum and therefore make possible to fabricate quantum de- 
vices with energy level spacing much greater than the temperature smearing kT at work 
temperature T (see, e.g., Ref. HJ). Moreover, dimension and shape of a quantum dot affect 
considerably the most important characteristics of the corresponding devices: relaxation and 
recombination time, Auger recombination coefficient etc, thus a possibility arises to control 
such characteristics in manufacturing the devices EIEHH Another way to control the properties 
of a quantum dot is instilling an impurity into the dot. Therefore, the investigation of spec- 
tral properties of a quantum dot with impurities as well as the dependence of the spectrum 
on the geometric parameters of the dot and physical characteristics of the impurity is an 
important problem of nano- and mesoscopic physics (see, e.g., in Refs. El El Eland references 
therein). The case of a hydrogen- like impurity is one of the most extensively studied up to 
now; however, the spectral problem in this case has no exact solution. On the other hand, 
short-range impurities can be investigated in the framework of the point potential theory 
(also called the zero-range potential theory). An important peculiarity of the point potential 
method is that the spectral problem for a point perturbed Hamiltonian is explicitly soluble 
as soon as the Green function for the unperturbed operator is known in an explicit formP^ 
For modelling the geometric confinement of a quantum dot, quadratic (in other words, 
parabolic) potentials are successfully usecP^ (see also examples of applications in Refs. 
13 ITTi -IT5|). The reason is that the self-consistent solution to the corresponding system of 
the Poisson and Schrodinger equations leads to the confinement potential having the form 
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of a truncated parabolic potential.^ Moreover, the Green function of the corresponding 
Hamiltonian H°, 

H° = -— A + ^-r\ (1) 

can be explicitly calculated^ E^lEO (here Q is the frequency of the oscillator, // denotes over 
the paper the mass of the considered charged particle). This makes possible to perform 
an exhaustive spectral analysis of the perturbation of H° by a point potential of arbitrary 
position q and strength a (we denote this perturbation by H a (q)) and to analyze the behavior 
of the eigenvalues of H a (q) as functions of q and a. This analysis is the main goal of the 
paper. Note that a quite particular case of the point perturbation of H° at q = (without 
obtaining any explicit form for the Green function) has been considered in Ref. 1213 Point 
potential for modelling an impurity in a spherically symmetric quantum dot has been studied 
in the series of papers using the Green function representation by means of the Laplace 
transform of the propagator kernel, but this approach allows to analyze (with numerical 
methods) the lowest impurity level onlyP^^^ 

It should be noted that point perturbations of the one-dimensional harmonic oscillators 
have been studied in detail earlier. This study was started in Ref. where the spectral 
properties of the point perturbed harmonic oscillator have been considered in the context of 
the one- dimensional models for the toponium physics and the Bose-Einstein condensation.^ 
A strict mathematical justification of results from Ref. l2*Tlwas done in Refs. I221I2U see also 
in Ref. 1213 Undoubtedly, our approach using the three-dimensional harmonic oscillator is 
more adequate for the analyzing the spectral properties of three-dimensional systems, in par- 
ticular, the toponium. It should be noted also that the one-dimensional harmonic oscillator 
perturbed by a point potential with varying position and strength has been investigated in 
Refs. [2H1 123 A series of phenomena of low- dimensional condensed matter physics can be ana- 
lyzed by means of the Hamiltonian of the perturbed oscillator: impurity in a one-dimensional 
quantum well, one-dimensional channel in a two-dimensional heterostructure subjected to a 
perpendicular uniform magnetic field etc., see the bibliography in the cited papers for details. 
However, the analysis given in Refs. |2Hl EH is based on the properties of one-dimensional 
second-order differential operators and is not extended to the three-dimensional case. 

The paper is organized as follows. Preliminary results are collected in Section |HI In 
Section Hill we consider point perturbations of the operator 

H° = ~A + V (2) 
2/i 

with an infinitely growing potential V. It turns out that the operator H a (q) can be defined 
and investigated for the more generic case when H° is defined by Equation (J2J). In Section ITVl 
some important properties of H a (q) are established. In particular, a complete description 
of the spectrum and eigenfunctions of H a (q) is given in Theorem 1. As a consequence 
of this theorem we get the falling of the considered particle on the attractive center as the 
potential strength a tends to — oo; for a very particular case of the one-dimensional harmonic 
oscillator perturbed at the potential minimum this phenomenon was observed in Ref. El 
In Section IS we define at fixed a a family of continuous functions such that the values of 
these functions at the point q form the complete family of the eigenvalues of H a (q). Some 
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elementary properties of these functions are established in Theorem 2. The main results of 
the paper are contained in Section IVIl where the point perturbations of the Hamiltonian 
of the harmonic oscillator are studied; the case of the isotropic harmonic oscillator {1} is 
considered in detail. These results are based on an explicit form of the Green function for 
the operator (JJJ. The detailed analysis of the dependence of the point levels on the position 
q and on the strength a is given in Theorem 3. In particular, if q ^ 0, then the point levels 
never coincide with the eigenvalues of the unperturbed operator H°. Therefore we have 
here no accidental degeneracy of the levels, which is a peculiarity of the one-dimensional 
model for the toponiump2Eia Hence, this degeneracy is an artefact of the one-dimensional 
model. Another interesting results is the asymptotic expression for the bound state of 
H a {q) (Equations (|3^jl. (JMJ)). These equations show that at least for the isotropic harmonic 
oscillator its potential (i.e., the frequency Q) can be recovered from the dependence of the 
ground state of the point perturbation on the support of the perturbation. Moreover, we 
argue that the form of the parabolic potential V may be recovered from the behavior of the 
excited energy for the ground state. Our conjecture is that this property is true for a more 
general form of the potential V. In this connection it is of interest to note that the study of 
the excited energy is one of the main problems of the quantum dot physicsPThe methods of 
Section IVTl allow us to analyze rigorously the phenomenon of so-called "positional disorder" 
in quantum dots (including nonisotropic ones). The relation of the degeneracy properties of 
the eigenvalues of -ffa(q) at q = to the symmetry properties of the unperturbed operator 
H° in the phase space is briefly discussed in the conclusion of Section IVIl In particular, the 
appearance of states with nonzero dipole momentum is noted. 



II. Preliminaries 

Here we present for the convenience of readers some basic properties of point perturbations 
of Schrddinger operators in L 2 (R 3 ) (see, e.g., Refs. 123 123 EBB ED for details). We will 
consider only Schrddinger operators H° of the form (J2J), where the potential V is subordinated 
to the conditions 

(PI) V E L p loc (R 3 ) for some p > 3; 

(P2) V- = min(V, 0) e L 2 (R 3 ) + L°°(R 3 ). 

Conditions (PI), (P2) are weaker than commonly used in applications conditions V G 
L^ C (R 3 ) and V > c with eel but making use of (PI), (P2) requires no change in proving 
of main results below. It is well known that under these conditions H° is semibounded from 
below and essentially self-adjoint on C£°(R 3 ) (see in Ref. 13*21 Theorem X.28). Further we 
put, as a rule, K = 1, /i = 1/2 and denote the obtained operator —A + V by H°. For the 
domain V(H°) of H° we have C^°(1R 3 ) C V(H°) C C(1R 3 ). This inclusion implies that the 
Green function G°(x, y; () for H° (the integral kernel of the resolvent R{() = (H° — C) _1 ) is 
a Carleman operator, this means that 

J |G°(x,y;C)| 2 rfy < +oo for a.e. xGK 3 . (3) 
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Moreover, according to Theorem B.7.2 from Ref. EHJ for every fixed (, ( G C \ spec {H°), 
the function G° obeys the following properties: 

(Gl) For every ( G spec(H°) the function G°(x, y; () is continuous in the domain {(x, y) G 
l 3 xl 3 : x ^ y}; 



(G2) |G (x,y;C)|<c 2 (C)|x-y| 



-l. 



(G3) if\x — y\ > d > 0, then |G°(x, y;()| < c 3 (d, 5, () exp(-5|x - y|) for some 5 > 0. 
Moreover, if Re ( < S = inf spec(H°), then arbitrary 5 with 5 2 /2 < £ — Re£ is 
suitable for this estimate. 

From (Gl) we get, in particular, that (jH)) is valid for every x G M 3 . 

The crucial role in the point potential theory is played by the regularized Green function 

G° cg (x, y; C) = G°(x, y; C) - ^ j^T ■ ( 4 ) 



In the particular cases, e.g. if V G C°°(R ), it is known that at fixed ( this function has a 
continuous extension on the whole space M 3 x M 3 (see, e.g., Ref. EHlor Theorem III. 5.1 in 
Ref. Elj). We need this property in the general situation and prove it under conditions (PI), 
(P2). 

It is sufficient to prove that G° eg (x, y; () is continuous with respect to (x, y) for some 
C = E Q < 0. Indeed, then for every ( G C \ spec (H°) 

C 

G°(x,y;C) = J ^G°(x,y;A)dA + G (x,y;£ ), 

where the path of integration lies in the resolvent set C \ spec (if ). The function 
(<9G°/<9A)(x, y; A) is jointly continuous with respect to (x, y) since it coincides with the 
integral kernel of (H° — A) -2 and this kernel is continuous according to Theorem B.7.1 from 
Ref. EU 

It is easy to see that V can be represented in the form V = V± + W, where V\ G C°°(1R 3 ) 
and obeys the property (P2) and W G L P (M 3 ) n L^M 3 ). Denote H 1 = -A + V u Si = 
inf spec (H l ) and by G 1 the Green function of H l . Fix E Q , E Q < min(S, S 1 ), and introduce 
the function F(x, y, z) = G°(x, z; E )W(z)G 1 (z, y; E ). Using properties (G2), (G3) and 
the estimate 

J |x-yr 

|y-a|<r 



where < v < 3, r > 0, a, x G M 3 , it is easy to prove that F(x, y, •) G L 1 (M 3 ) for all 
x, y G M 3 . In virtue of the Lippmann-Schwinger relation 

G°(x,y;£ ) = G' 1 (x,y; J E; ) + f G°(x, z; E )W(z)G\z, y; E ) dz 
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and the continuity of the regularized Green function for H 1 , it remains to prove that the 
function 

J ( x >y) = J F ( x ,y, z ) dz 

is continuous on R 3 x R 3 . Moreover, (Gl) shows that it remains to prove the continuity 
of / at points of the form (xo,xo). To do this fix e > and find rj > such that the 
relations |x — x | < 77, |y — x | < r\ imply |/(x, y) — I(x ,x )| < e. Introduce the sets 
Bi(.v) = { z : l z ~~ x o| < v}> B 2 {rj) = R 3 \ -81(77), and for a measurable set B Cl 3 denote 
7 B (x, y) = J B F(x, y, z) dz. Then 

l^(x,y) -/(x ,y )| < |/ Bl (7,)(x,y)| + \I Bl ( v )(^o,yo) \ + l^ 2 (??)( x >y) - ^(^(xo, Yo)| • 

If x,y,zG (77), then by (G2) 

|F(x,y,z)|</(z) (x-yl^lz-yr 1 

where / G L p , therefore relation (jSj) and the Cauchy-Schwartz inequality lead to the estimate 
|/ Bl (^)(x,y)| + |/ Bl (^(x ,yo)| < const?]. On the other hand, if x,y G B 1 (r)/2), z G B 2 {fj), 
then we have from (G3): |F(x, y,z)| < g(z) exp(— S\z\), where 5 > and g G L p . Thus 
by (Gl) and the Lebesgue majorization theorem, J B2 ( r? )(x, y) is a continuous function on 
-81(77/2) x B 1 (i]/2), and the proof of continuity of G° eg is completed. 

Let q G M 3 , then the restriction of H° to the domain {/ G V(H°) : f(q) = 0} is 
a closed symmetric operator 5* with the deficiency indices (1,1). By definition, the point 
perturbation of H° , supported on q is a self-adjoint extension of 5* different from H°. All the 
point perturbations of H° supported on a given q G M 3 form a one-parameter family H a (q), 
a G M, of self-adjoin operators such that the Green function G a of H a (q) is given by the 
formula 

G a (x, y; C) = G°(x, y; C) - [Q(C; q) - G°(x, q; C)G°(q, y; C) , (6) 

which is a consequence of the Krein resolvent formula. Here Q((; q) = G° cg (q, q; C) is the 
so-called Krein Q-function. The operator H° corresponds formally to a = 00; moreover, H° 
is the Friedrichs extension of S. 

The extension parameter a has an important physical meaning, namely, H a can be 
treated as the Hamiltonian H° perturbed by a zero-range potential, in this case a is the 
strength of this potentialP^"^ In place of the strength a, it is more convenient to use for 
applications so-called "scattering length" £ s , £ s = l/(47ra) (see in Refs. |H1 ESI HUD again) . 
More precisely, 

t =-^- 

s 2nh 2 a ' 

and we see that l s has actually the dimension of the length. 

Note that according to the general results of the Krein self-adjoint extension theory, 
the function ( 1— > Q(£; q) is analytic in the domain C \ spec(H°) for each q G M 3 and 
dQ(E; q) / dE > if E G R \ spec {H°) W Remark that Q((; q) can be continuously extended 
to some points of spec(if°). Further we assume that Q((] q) is continuously extended to all 
regular points. 
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It is easy to prove that for every q G M 3 the mapping ( i— > G°(-,q;£) is an analytic 
function from the domain C \ spec (H°) to the Hilbert space L 2 (M 3 ). Denote G°(-, q; () by 
^q(C)j then we can rewrite (JBJ) in an operator form 

R a (0 = R°(C) - [Q(C; q) - a]" 1 |<? q (C)}<Sq(C)| , (7) 

where H«(C) = (# Q - C) _1 and R°(() = (H° - C) _1 . 

Note, that g q (C) is a nonzero function for every qGl 3 and £ G C \ spec Indeed, 
otherwise we have y(q) = for every G V(H°) that contradicts the inclusion C^°(IR 3 ) C 
V(H°). 

In conclusion we mention a possibility to approximate the zero-range perturbation by 
potentials with decreasing support. For V = the corresponding procedure is described in 
Ref. El (Theorem 1.2.5). We sketch here the proof for H° with potential V having properties 
(PI), (P2). 

Let W G L^ omp (R 3 ), in particular, W is a Rollnik function (see in Ref. IH21 Sec. X.2). 
Denote v = IW] 1 ^ 2 , u = v sign(V^), and let \(e) be a real-analytic function in a neigh- 
borhood of zero such that A(0) = 1. For e > consider the operator H £ = H £ (q) = 
H° + e- 2 A(e)W(£ _1 (x - q)). Then the resolvent R £ (() = (H £ - C) _1 (e > 0) has the form 

R £ (0 = R°(0-eX(e)A £ [l + B £ ]- 1 C £ , 

where A 6 , B £ , C £ are integral operators with the kernels A £ (x, y; £) = G°(x, ey + q; Qv(y), 
C £ (x,y;C) = G (£x + q,y;C)M(x), 5 £ (x, y; C) = eA( £ )G (ex + q, ey + q; ()u(x)v(y). Define 
A and C° putting e = in the formulas above, and define B° by the integral kernel 
-B°(x, y) = (47r|y — x|) _1 n(x)t>(y). All the operators A £ , B £ and C £ (e > 0) belong to 
the Hilbert-Schmidt class and A £ — > A , B £ — > Z?°, C £ — > C° with respect to the Hilbert- 
Schmidt norm as e — > +0. Moreover, using (jlj) we can prove that with respect to this 
norm 

B £ = B° + e(\'(0)B° + Q(C,ci)\u)(v\) + o(e) . 

Hence, the arguments using for the proof of Theorem 1.2.5 from Ref. [HI give the following 
result. 

Theorem A. 

(1) Let (v\ip) = for all L 2 -solutions ip of the equation B°ip = —if (in particular, let —1 be 
not an eigenvalue of B ). Then H £ (q) — > H° in the norm-resolvent sense as e — > +0. 

(2) Let —1 be a simple eigenvalue of B° and if be a corresponding eigenf unction normalized 
by the condition (<f\if) = — 1, where <f> = <fsign(V). If (v\if) ^ 0, then lim H £ (q) = 

H a (q) in the norm-resolvent sense, where a = —\'(0)\(v\if)\~ 2 . 

(3) Let —1 be a multiple eigenvalue of B° with eigenf unctions (fi, . . . , ip n normalized by the 
conditions (ifj\ip k ) = —Sjk (<fj = fj sign(l^)). If (v\(fj) ^ for some j and A'(0) ^ ; 
then lim H £ (q) = H a (q) in the norm-resolvent sense, where 



a = -A'(0) 



.i=i 
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III. Point perturbation in the case of unbounded po- 
tential V 

Starting with this section we suppose additionally that 
(P3) lim V(r) = +00. 

|r|^oo 

In this case -R°(C) is a compact operator for all £ G C\spec(if ) (the Strichartz theorem; see, 
e.g., in Ref. EHl Theorem XIII. 69). Therefore spec(if°) consists of an unbounded sequence 
Ao < Ai < . . . < \ n < . . . of eigenvalues with finite multiplicity k n . Consequently, Q((; q) is 
a meromorphic function of (. We are going to find the poles of this function. 

Denote by L n the eigenspace associated with A n , and choose in L n an orthonormal basis 
F„ jfc (r), k = 1, . . . , k n . For every q G M 3 we denote 

a(q) = {X n G spec(#°) : 3f G L n s.t. /(q) + 0} 



Lemma 1. The set of all poles of the function ( 1— > Q((; q) coincides with <r(q). 

Proof. Since (dG° /d()(x, y; () is the integral kernel for the operator (H° — ()~ 2 , we have 
according to the Mercer theorem 




where the series converges locally uniformly on I 3 x R 3 x (C \ spec(i/ )). Therefore, 

r\ OO k n 

^Q(C; q) = E E ( a « " 0- 2 l^, fe (q)| 2 , (8) 

^ n=0 fe=l 

and the series converges locally uniformly on (C \ spec(H )) x R 3 . The lemma follows from 
(jHJ) immediately. I 

Lemma 2. For each q G M 3 the set <r(q) is infinite. If V is bounded from below, then 
Ao e tr(q). 

Proof. Consider the space of continuous functions C(1R 3 ) with the topology of compact 
convergence. Due to the closed graph theorem and the relation V(H°) C C(1R 3 ), the operator 
R°(— 1) : L 2 (R 3 ) — » C(1R 3 ) is continuous. Therefore for every / G V(H°) the Fourier 
expansion for / with respect to the basis {F n ^) n ^ converges locally uniformly. Assume that 
the set cr(q) is finite; let = maxjn : A„ G o"(q)} and P be the orthogonal projection 
of L 2 (M 3 ) on the subspace M = Lq + . . . + Ln- Then for every tp G T>(H°) the conditions 
</2(q) = and (P(p)(q) = are equivalent. Since M is finite dimensional, there is h G M 
such that for every (p G M the conditions <£>(q) = and (h \ (p) — are also equivalent. Using 
the inclusion (^(M 3 ) C V(H°) we see that there is a function h G L 2 (IR 3 ) such that for 
every ip G (^^(IR 3 ) the conditions ^(q) — and (h\<p) = are equivalent. Obviously, this 
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is impossible, hence cr(q) is infinite. If V is bounded from below, then by Theorem XIII. 48 
from Ref. EH the eigenfunctions of H° corresponding to the ground state Ao have no zeros 
therefore A G a(q). I 

Another property of the function £ i— » Q(C', q) we need further follows. 
Lemma 3. The function Q((; q) tends to — oo as ( — > —oo, ( G R. 

Proof. Since if is the Friedrichs extension of S 1 , the statement follows from Proposition 4 
of Ref. EHl I 

IV. Spectral properties of H a at fixed position of the 
point perturbation 

Here we describe the spectrum of H a (q) for a fixed q G R 3 . Further, if it does not lead to a 
misunderstanding, we omit q from the notations. 

Since H a is a rank one perturbation of H°, the spectrum of H Q is discrete. Moreover, an 
eigenvalue A n of H° of the multiplicity k n is an eigenvalue of H Q of the multiplicity k n — 1, 
fc„ or fc n + 1 (if k n = 1, the first case means, of course, that A n does not belong to spec(H a )). 
For A ^ spec(H°) we see from (J7J) that A is an eigenvalue of H a if and only if ( = A is a 
solution to the equation 

Q(C;q)-a = 0. (9) 

Denote by (e n )neN = (£n(q))neN the strictly increasing sequence of all the poles of q). 
Since (dQ/dE)(E; q) > for £J G R\spec(if°), the equation Q has exactly one solution on 
each interval (— oo,£ )j (^o^i); • • • Denote such solutions, which do not belong to spec(H°), 
by S ,Si, . . ., where £ < £\ < ■ ■ -. The following theorem completely describes the eigen- 
values and the eigenfunctions of H a (q). 

Theorem 1. Let q G R 3 be fixed. The spectrum of H a = H a (q) is discrete and consists of 
four nonintersecting parts o~i, a 2 , <7 3 , cr 4 described as follows. 

(1) o\ is the set of all solutions S n to the equation (0), which do not belong to spec(H°). 
The multiplicity of £ n in the spectrum of H° is equal to 1. 

(2) o"2 is the set of all A n G cr(q) that are multiple eigenvalues of H . The multiplicity of 
the eigenvalue X n G cr 2 in the spectrum of H a is equal to k n — 1 . 

(3) (T3 consists of all X n , A n G spec (H°) \ cr(q), that are not solutions of (0). The multi- 
plicity of the eigenvalue A n in spec(H a ) is equal to k n . 

(4) cr 4 consists of all X n , A n G spec (if ) \ cr(q), such that X n is a solution of (0). The 
multiplicity of the eigenvalue X n in spec(H a ) is equal to k n + 1. 

The corresponding eigensub spaces are described as follows. 
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(1) The subspace spanned by the normalized eigenf unction 

i 

2 

#q(£n) • 

(2) The orthogonal complement in L n of the function 

fen 

* n (x) = ^F„ ifc (q)F n , fc (x) , 
fc=i 

or, equivalently, the subspace of L n of the form {f £ L n : /(q) = 0}. 

(3) T/ie subspace L n . 

(4) T/ie direct sum of L n and the space spanned by the function o q (A n ), which is orthogonal 
to L n . 

Proof. The proof is based on direct calculations with the help of following statements: 

(A) The orthoprojector P(Eq) on the eigenspace of a self-adjoint operator T corresponding 
to an isolated eigenvalue E has the form 

P(E ) = -Res[(T - Cr X ;C = E }. 

(B) Suppose P11P2 and P\ + CP2, where c G C, are orthoprojector s in a Hilbert space and 
P2 7^ 0, then c equals 0, 1 or —1. 

The first statement is well-known; we omit the easy proof of the second one. Denote by 

MO, 

^(C) = [Q(C;q)-«]- 1 |a q (C))(^(C)l, 

the second term in the representation (JJJ) of the resolvent. Further, denote for E G R 

P a (Ea) = -Res[Ra(();( = E \, 

P°(^) = -Res[i2°(C);C = ^o], 
T(£ ) = Res[A(C);C = £o]; 

therefore according to (JJJ) 

P a (E ) = P°(E )+T(E ). 

Start with the proof of the first assertion of Theorem. It is obvious that o\ C spec (H a ). 
Let £ n G (T\, then in a vicinity of £ n we have the following expansion 

Q(C; q) - a = ^-Q(S n , q)(C - S n ) + 0(( - E n f . (10) 

Therefore 

|s q (£ n ))(<7 q (5 n )|. (11) 



T{E n ) 



d 

— Q(£ n ;q) 
9C 
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Since obviously P°(£ n ) = 0, we have P a (£ n ) = T(£ n ) and the normalized eigenfunction 
corresponding to £ n is 

i 

2 g^Sn) ■ (12) 

Now consider an eigenvalue A n of if . In this case -P a (A n ) = P°(A n ) + T(A n ). According 
to (jHj), in a neighborhood W of X n we have the following representation 

( ?q (C) = ^„(-;q)(A„-C)- 1 + /(C), 
where / is analytic function in W with values in L 2 (R 3 ) and 

^„(x; q) = ^F nifc (q)F nifc (x) . 
k=l 

Consider the following three cases: (a) A n G cr(q); (b) A n ^ cr(q) and Q(A n ;q) - ft^O; 
(c) A n ^ cr(q) and Q(A n ; q) - a = 0. 

Let us start with the case (a). Since A n is a pole of Q(-;q), we have ^ n (-;q) ^ and 
therefore T = cP, where P is the orthoprojector on the one-dimensional space spanned by 
\l/ n (-;q). Since \I/ n (-;q) G L n , in virtue of statement (B) c = —1, and the assertion (2) of 
Theorem is proven. 

In the case (b) according to Lemma 1, F n ^(q) — for all k = 1, . . . , k n ] hence \l/ n (-, q) = 
and T(A n ) = 0. This implies assertion (3) of Theorem. 

Finally in the case (c) we can use (JTUJ) — (H2J) with ( = \ n instead of ( = £ n , and obtain 

T(A n ) = |$ n )($„|, 

according to (B), this get the statement (4) of Theorem. I 

For n G N denote by A n the set of all a G R such that the solution £ n = £ n (a) of 
Equation (JHJ) does not belong to the spectrum of H°. Lemma 2 shows that R \ A n is finite, 
moreover, if V bounded from below, then A = R. 

For all q G R 3 we will denote £_i(q) = A_i = — oo. Using Lemmas 1 and 3 we get 
immediately the following proposition. 

Proposition 1. For each n G N the function a t— > £ n {ct) strictly increases on A n . Moreover, 

lim £ n (a) = e n , lim £ n (a) = e n -\ .1 

Remark. For n = 0we have an interesting phenomenon of falling the considered particle on 
the point q (the falling on the attractive center; cf. Ref. |^]for the case of a one-dimensional 
oscillator). Indeed, using estimate (b') from Theorem B.7.1 of Ref. we obtain without 
any difficulty |$o( x )| 2 ~~ > 5(x — q) in an appropriate space of distributions as a — >• — oo (and 
therefore £q — > — oo). According to the standard interpretation of quantum mechanics, this 
relation means that the probability to find the particle in a domain not containing the point 
q tends to zero as £q tends to — oo. 



-^r(£n;q) 
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V. Dependence of the spectrum of H a (q) on q 

Here we are going to analyze the dependence of the eigenvalues of H a on q. It is clear 
that S n (q) are continuous branches of the multi- valued function defined by equation (JHJ). 
This branches can intersect at values A n where a monodromy arises. To get a univalent 
enumeration of these branches, we modify the parameterization of the eigenvalues of H a 
given by Theorem 1 (the enumeration of the numbers S n (q) depends on the enumeration of 
poles e n G spec(H°), which in its turn depends obviously on q). For n = —1, 0, . . . consider 
the sets X n defined as follows: X_i = R 3 , and 

X„ = {q G R 3 : 3/ G L n s.t. /(q) ^ 0} = {q G I 3 : A n G a(q)} , 

for n > 0. For all n G N the set R 3 \ X n is nowhere dense in IR 3 (see in Ref. EH1 Theorem 
XIII. 63). According to Lemma 1, for n > 0, the set X n coincides with the set of all q6 R 3 
such that A n is a pole of the function q). Since we do not suppose the potential V is 
smooth, the function Q((;q) on the set (A n _i,A n ) x (X n _i fl X n ), n > 0, is not, generally 
speaking, smooth. Nevertheless, it is monotone and real analytic with respect to the first 
argument ( and continuous with respect to the second argument q. In this case the following 
simple variant of Implicit Function Theorem is applicable (see in Ref. HHfor the proof): 

Let J be an open nonempty interval of the real line R ; X be a topological space, and F : 
J x X — > R be a separately continuous function such that each partial function t i— > F(t,x), 
x G X , is strictly monotone. Suppose that F(t ,xo) = for some (t ,xo) G J x X. Then 
there are an open neighborhood U of the point Xq in X and a continuous function f : U — > J 
such that (1) F(f(x),x) = for all x G U; (2) if U' is another neighborhood of Xq, and 
g : U' — ► J is a function with the property: F(g(x),x) = for all x G U' , then U' C U, and 
f\u> = 9- 

According to this version of Implicit Function Theorem, for any q G X n -\ fl X n there 
exists a unique solution E n (q) to equation that belongs to (A n _i, A„) and q i— ► E n (q) is 
a continuous function in X n _\ fl X n . 

Proposition 2. Every function E n (q), n = 0,1,..., has a continuous extension to the 
whole space R 3 . 

Proof. Fix n = 0, 1, . . ., and let a point q, q G R 3 \ (X n _i fl X n ), be given. Choose 
a sequence (qfc)feeN from X n _x fl X n which tends to q. First we note that the sequence 
(E n (qk))keN is bounded in R. It is trivial for n > 0. If n — 0, the sequence is bounded 
from above. We prove that it is bounded from below as well. Otherwise Eq{c^) — » — oo for 
some subsequence (q^J. Since Q(E; q) — > — oo as E — > — oo, there exists A < A such that 
Q(A; q) < a. Then there exists iV G N such that Q(A; q^) < a and i^oCoO < A if I > N. 
Therefore, for k > N we have 

Q(E (q kl ); q*,) - « < q fei ) - a < , 

and we get a contradiction with the definition of £/o(qfcJ- 

By Bolzano- Weierstrass we can extract a subsequence (q^) from the sequence (q^) such 
that the subsequence (-E^q^)) has a limit, which we denote by E' . To prove that the 
sequence (E n (q k )) tends to E' and E' is independent of the choice of a sequence (q^) tending 
to q we need the following lemma concerning properties of E '. 
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Lemma 5. The limit E' has the properties: 

(1) E' is not a pole of the function ( i— > Q((; q); 

(2) if \ n -i < E' < \ n , then E' is a unique solution of Equation Q in the interval 
(A n -i, \ n ); 

(3) if E' = A„_i, then lim [Q(£; q) - a] > 0; 

E — >E' 

(4) = \ n , then lim [Q(£7; q) - a) < 0. 

E — >E f 

Proof of the lemma. 

(1) First consider the case n > 0. The function Q„(C; q) = [Q(Ci q) — a ](C — A n _i)(C — 
A„) is continuous on the interval (A n _ 2 , A n+ i) x IR 3 . Since Q n (E n (q kl ); q kl ) = 0, passing 
to the limit I — ► oo we get Q n (E';q) = 0. Suppose ( = E' is a pole of Q(C;q), then 
Q n (E';q) = Res[Q(C;q);C = E'} 0> an d we a contradiction. For n = 0, we consider 
Qo((] q) = [Q(Ci Q) — a ](C ~ Ao), and get the same result. 

(2) It is sufficient to pass to the limit / — > oo in the identity Q(E n (q kl ); q/cj = 0. 

(3) In virtue of statement (1) of the lemma, the function ( i— > Q((;q) is continuous in 
a neighborhood of E', and therefore there exists a limit lim [Q((; q) — a] = L. Assume 

that L < 0, then Q(-E,q) — a < for some i? G (A„_i,A„). Choose some m such that 
E n (q km ) < E. Since Q((; q) increases on the interval (A„_i,A„) as the function of (, we 
obtain a contradiction: 

= Q{E n (q km )- q km ) -a< Q(E; q k J - a < . 

Statement (4) can be proven similarly to (3). I 

Let us return to the proof of the proposition. We prove that if a sequence (p k ) k e?q from 
Xfc_i fl X k converges to the point q, then E n (p k ) — > E'. 

Suppose E n (p k ) does not converge to E', then there exists a subsequence (p^J such that 
E n (p kl ) — > E*, E* 7^ E' . Assume E* < E' . Taking into account item (2) of Lemma 5 we get 
E* = A n _i or E' = X n . In both the cases we have 

lim [Q(C; q) - a] > and lim [Q(C; q) - a] < 0. 

£ — £ — >E' 

Take some real numbers E\ and -E2 such that E* < E\ < E2 < E'. Then by the strict 
monotonicity of £ 1— ► Q(C; q) we have 

< Q(Ei; q) - a < Q(£ 2 , q) - or < . 

This is a contradiction. I 

The following theorem is the main result of this section. 

Theorem 2. For each fixed a G R £/iere is a sequence (E n (q)) ne ^ of continuous functions 
0/qGl 3 with the following properties: 
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(1) A n _i < E n (q) < X n for all n G N. 

(2) For each qe I 3 the set consisting of all E n (q) and all the numbers X n with multiplic- 
ities k n > 1 form the complete collection of the eigenvalues of the operator H a (q) . 

(3) If A n _i < E n (q) < X n , then E n (q) is a unique solution of the equation (0) on the 
interval (A n _i, A n ). 

(4) If ( = X n is a pole of the function ( i— ► Q((] q), then E n -\(q) < X n < E n (q). 

(5) If ( = X n is not a pole of the function ( i— > <5(C;q) ; then we have the following 
assertions: 

(a) ifQ(X n ; q) - a < 0, then E n (q) = X n < E n+1 (q); 

(b) ifQ{X n ; q) - a > 0, then E n (q) < X n = E n+1 {q); 

(c) ifQ(X n ; q) - a = 0, then E n (q) = X n = E n+1 (q). 

Proof. Consider the functions E n (q) given by Proposition 2. Then (1) is obvious by 
definition of ^(q). Assertion (2) follows from Theorem 1. Assertions (3) and (4) were 
proven in Lemma 5. It remains to prove (5). 

Let A n be not a pole of ( i— > Q((; q). Suppose Q(X n ; q)— a < 0. For any positive integer m 
we choose a number E' m such that X n — 1/m < E' m < A n ; then Q(E' m ; q) — a < 0. Further, we 
choose points q m G M 3 such that A n _i and A n are not poles of the function ( i— > Q(C', q m ) (that 
is q m G A m _inA m ), and such that |q— q m | < 1/m and Q(E' m ; q m )— a < 0. Then ( = E n (q m ) 
is a solution of the equation Q(C', q m ) — a = lying in the interval (A n _i, A n ). Since Q((; q m ) 
is a strictly monotone function of ( on this interval, the inequalities E' m < E n (q m ) < X n take 
place for all m. Thus E n (q m ) — > X n and q m — > q as m — > oo; therefore A n = £7 n (q) by the 
definition of the function £J n (q). According to Lemma 5, Q(A n ; q) — a > 0, if X n = E n+ i(q); 
therefore A n < E n+ i(q). Hence, item (5a) is proved. The proofs of items (5b) and (5c) are 
similar. I 

Theorem 2 gives a useful description of the spectrum of H a . Namely, denote by M the set 
{m G N : k m > 1} and together with the functions E n (q) introduce a sequence of constant 
functions A^(q) = A m , where m G M, k — 1, . . . , k m — 1. Then i? n (q) < A^(q) < E ri+1 (q) 
for all n G M, fc = 1, . . . , k n — 1, and for any fixed q G M 3 the union of the sequences 
(E n (q)) nm and (A^ ) (q)) me A/,fc=i,...,A :m -i forms the complete set of the eigenvalues of H a (q) 

oo 

multiplicity counting. If q G O X n , then every -E n (q) is distinct from the numbers A^(q). 

n=0 

oo 

Since E \ P| A n is the set of the first Baire category, for a generic q the point perturbation 

n=0 

levels -E n (q) are distinct from the levels of the unperturbed operator H°. 
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VI. Point perturbations of the harmonic oscillator 



Here we apply the results of the previous sections to the Hamiltonian (J2J) with the potential 

V(r) = ^ + ^ + ^, (13) 

where Qj (j = x, y, z) are the frequencies of the oscillator. The function V can be considered 
as a confinement potential of a quantum well in M 3 with the characteristic sizes 



Lj -^2^V i- x >v> z > 

(numbers V%Lj are called also length parameters of the oscillatoJ^) . Therefore the operator 
with potential (jlHj) can be used as the Hamiltonian of a (generally speaking, asymmetric) 
quantum dot It is convenient to pass to dimensionless coordinates x = r/L, where L = 
^jL x L y L z . In the coordinates x = (xi,X2 ) x%) the operator H° takes the form H° = h~QH°, 
where 

zjO a i 1 / 2 2 i 2 2, 2 2\ 

H = -A + - (uj 1 x 1 + uj 2 x 2 + uj 3 x 3 ) , 

f2 = ^/Vt x VL y Vt z , loi = — , u>2 = , ^3 = -^y 
(hence, wia^a^ = 1). 

Further we discuss the properties of if . The spectrum of this operator consists of the 
eigenvalues 

Anin 2 n 3 = + V 2 ) + ^(™2 + 1/2) + ^ 3 (™3 + V 2 ) , 

where ni, n2,n 3 G N. The corresponding normalized eigenfunctions are 



(x) = ^1(^1)^2(^2)^3(^3) 



where 



^fe)=(g) 1/4 (2^.!)-^exp 




is the oscillator function (H n (x) is the Hermite polynomial of degree n). 

If the frequencies oo%, 002,003 are independent over the ring Z (this is the generic case), 
then the spectrum of H° is simple; therefore the multiplicity of the eigenvalues of H a (q) 
does not exceed 2 and the part 02 of the spectrum spec (H a (q)) is always empty. On the 
other hand, since H n (0) = if and only if n is odd, X ni ,n 2 ,n 3 £ spec(if Q (0)) if and only if one 
of the numbers rij (j = 1, 2, 3) is odd; hence, spec(H°) \ cr(0) is always infinite. In addition, 
for all n > the set M. 3 \ X n is infinite. 

In general case, there are no explicit expressions for the Green functions of the harmonic 
oscillator in terms of commonly used elementary or special functions. Nevertheless, in a 
number of cases, the representation of the Green function G°(x, y; E) as the Laplace trans- 
form of the heat kernel K(~x, y; t) for H° is very useful to investigate some properties of the 
Krein Q-function. The heat kernel for H° has the form (see, e.g., in Ref. 14*2")) : 
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Using the heat kernel K* for the free Hamiltonian H* = — A, 

^(x,y;t) = (4vrt)- 3 / 2 exp ("^^) 



and the Q-function for H 



Q f (C) 



47T 

we get immediately from the formula 



G(x,y;£) = J e m K(x,y;t)dt, 
o 

that for Re ( < (u>i + u> 2 + uj 3 )/2 the following representation of the Q-function for H° takes 
place: 



Q(C;q) 



-C 



47T 




It is clear from (JHJ) that (dQ/dqj)(E; q) < for ^ > 0, if £ < A = (wi + u 2 + w 3 )/2. 
Since > for E E R \ spec (H°), © implies that dEo/dqj > 0. In particular, the 

depth of the lowest impurity level A — E (q) decreases if |q| increases in such a way that 
the inner product a • q remains positive for each vector a with positive coordinates. In the 
spherically symmetric case U\ = oj 2 = ^3, we have dQ/dq < 0, where q — |q| > 0, and the 
depth decreases with increasing of q. This phenomenon was discovered numerically for a 
spherically symmetric quantum dot in Ref. El and called positional disorder. We see that 
the positional disorder is common to each parabolic quantum dot, not only to the spherically 
symmetric one. The similar result is valid in the two-dimensional case, i.e. for the case of 
impurities in a quantum well (see numerical results in Ref. IHj) . Our arguments are valid 
in the two-dimensional case also, thus we have a strict proof for the positional disorder in a 
two-dimensional quantum well. 

The more detailed analysis is possible in the case of the isotropic oscillator. Q x = Q y = Q z 
(= Q), i.e. in the case of a spherically symmetric quantum dot. In this case uj\ = uo 2 = 003 = 1 
and the spectrum of H° consists of the eigenvalues 

3 

A n = n + - , neN, 

where X n has the multiplicity k n = (n + l)(n + 2)/2. In this case there are natural units 
of length (namely, L) and of energy (hQ). Therefore the following very important scaling 
properties takes place. Denote by Q(Ci o) t ne Krein Q-function for the operator H° keeping 
the notation Q(C; q) for the Q-function of H°. Then 
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Denote fi/(2nh 2 L) by a ; obviously, a is strength of the point potential corresponding to 
the scattering length L. Then Equation © takes the form 



or, equivalently, 



Equation (|15p shows that a change of the frequency Q does not change the numerical values 
of energy levels in the spectrum of H° if L is used as the unit of length, Ml as the unit of 
energy and a as the unit of point potential strength. 

In the case of isotropic oscillator, the set cr(q) has a simple description: 

Proposition 3. Let Q x = Q y = Vt z . Then cr(q) = {A 2 n '■ n £ N}, if q = 0, and er(q) = 
spec(-£T°) otherwise. 

Proof. Each A n is equal to A ni „ 2n3 , where n\ + n 2 + n% = n. If n is odd, then at least one 
of rij is odd, and ^ nin2n3 (0) = 0. Therefore A n ^ cr(0). On the other hand, if n is even, then 
^nOo(O) 7^ 0, and therefore A n G cr(0). 

Let now q ^ 0. First we remark that for all n G N the following assertion is valid: 
Lemma 6. If H n (x ) = 0, then H n+1 (x ) ^ 0. 

Proof of the lemma. For all n G N the following relation takes placed 

H' n+1 (x) = 2(n + l)H„{x). 

If H n (xo) = H n+ x(xo) = 0, then H' n (xo) = 0. Since y = H n (x) is a solution to the differential 
equation y" — 2xy' + 2ny = 0, we have H n [x) = for all x; but this is impossible. I 

Let us return to the proof of the proposition. Suppose that q ^ 0; without loss of 
generality we can assume g 2 ^ 0. Since H\(x) = only for x = 0, and H (x) ^ for all x, 
we have A , Ai G cr(q). Let n > 1. Suppose that $ n _i j i i0 (q) = 0, then according to Lemma 6, 
$n,o,o(q) + 0. I 

Using Proposition 3 we can give the complete description of the spectrum H a (q) in the 
case of an isotropic H°. Moreover, in this case the explicit form of the Green function 
(j°(x, x';£) is known, and therefore we can give the explicit form of the Krein Q-function 
and eigenfunction of H a (q). In particular, the equation for the point perturbation levels 
E n (q) can be obtained in an explicit form. The mentioned Green function has the formP^ 

'U(-(; 0U'{-C -V) + U'{-C, 0U(-(; -rj) 



^(x,y;C) = -^r(i- C 



, ^(-C; 0^(-C; -v) - ^(-C; 0^(-C; -v) 



X 



(16) 



where ^ = (|x + y| + |x — y|)/2, rj = (|x + y| — |x — y|)/2, U{y; z) is the parabolic cylinder 
functiorP^ (in the Whittaker notation U{v\ z) = D_ v _i/ 2 (z)), and U' denotes the derivative 
of U with respect to the second argument 

U'(C;y) = § y U(C;y). 
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Using ()16|) . we get the following expression for the Q-function: 

x [ (q 2 - O U(-(, q)U(-(, -q) + 4U'(-(, q)U\-^ -q) 

- - q {U'{-C q)U(-C, -q) - U(-C, q)U\-C, -q))], (17) 

where q = |q|. Due to the symmetry of the problem, the Q-function depends on q only, 
so we shall write often Q((;q) instead of Q{(;q)- Introducing the notation U((;y) = 
U((] y)U((; —y), we can rewrite (|T7j) in the sometimes more useful form 

Q(C; q) = -TT^rrr (\ - c) [ (g 2 - 4C) W(-C; g) - -W '(-C; g) - U "(-C; g) 

4(27r)2 \z / l q 

where the prime denotes the derivative with respect to the second argument as before. 
Passing to limit we get at q = 

1 I 4 2 

Q(C;0) = ~ T^v- < 19 > 

Vo7T -n / 1 S> 



It is interesting to compare (HH with the Krein Q-function Q (1) (C; 0) for the one-dimensional 
harmonic oscillator)^ 



qW(C;0) = 2-^ 



4 2 
4 2 



Curiously, in the case of the free Hamiltonian H° = —A, the Q-functions Qd for d = 1 and 
for <i = 3 are also related as follows: 

gr 1 (C) = -87rg 3 (C)- (20) 

Namely, for the free Hamiltonian Q^Q = (2 v ^()-\ Q 3 (() = -(47r) _1 >/ = £. For q ^ 
relation (|2"UJ) for Q-functions of the harmonic oscillators is violated. 

It is useful to consider the behavior of the function ( i— > Q((;q) near the singular points, 
i.e., near the poles and in a neighborhood of — oo. Using properties of the parabolic cylinder 
functions,^ we have 

Q(C; o) = - f -V - in 2 + 1 - hTT9^ + - x ^ 

(27r)2(2n)!! - A 2 „ 2 ^ fc(l + 2fc) 

as C — >■ A2n- If q 7^ 0, the coefficients for corresponding asymptotics are cumbrous enough, 
and we give the leading term only: 

exp (— g 2 /2) 
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+V2 (q- 1 - q) H n (g/W) H n+1 (q/V^j + H 2 n+X (g/v^)) (C - K)' 1 + 0(1) , (21) 



as C — > A n . 

For Re( — > — oo, we have 

Q(C,q) = - 



Arc 



(22) 



coincides with the Krein Q-function for 



It is important to note that the leading term in 
the free Hamiltonian —A. 

Now consider the properties of the function q i— > Q{(;q)- Since U(v;z) is an entire 
function of z, the function q \— > Q((; q) at ( spec(if °) can be extended to a real analytic 
even function on IR. (see (I18|0 . In particular, 



d_ 

dq 



Q(C;0) = o. 



As to the second derivative, we can obtain after some algebra 



dq 



2 Q(C;0) 



l 



8^' 



7T 



(4c 2 + i) l't 2 -sc u 2 



(23) 



For the fixed £ G K \ spec (if ), the asymptotics of Q at q — > oo is given by 



87T 



2C 1 + 2C 2 , _ 

g = — + 

q q 6 



(24) 



This follows from the asymptotics for U{C\ q) at q — > ooP^ 

a/2T 



W(C;?) 



r(| + C) 



where X = \/g 2 + 4£. 

Further the following formula will be also useful 



f«;e> 



rfl- 



— = — r — Fo ( ^ - C 

4v/2vrr(i-r l ° 



(25) 



Here and below we use the standard notations^ 



V(z) 



r^) 
rr*) 



The plot of the graphs for the function Q((; q) is shown on Fig. ^and Fig. El 

In the case of an isotropic oscillator, the functions E n (q) depend only on q and we will 
denote them by E n (q). Further properties of these functions (and, in particular, of the 
spectrum of H a (q)) for the isotropic case are given in Theorem 3 below, which is one of the 
main results of the article. 



Theorem 3. The following assertions take place. 
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Figure 1: Q as a function of £ for (a) g = 0, (b) q — 1/10, (c) g = 1, (d) q = 3. 
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(1) The functions E n (q), n G N are real- analytic. If a = and n > 0, then in a vicinity 
of zero, these functions are continuous branches of a two-valued analytic function. 

(2a) E (0) < A for each a, aeR. 

(2b) lfa>0, then E 2n+1 (0) = X 2n+1 and X 2n+1 < E 2n+2 (0) < A 2 „ +2 Vn G N. 

(2c) lfa<0, then A 2 „ < E 2n+1 (0) < A 2n+2 and E 2n+2 (0) = \ 2n+2 \/n G N. 

(2d) If a = 0, then E 2n+1 {0) = E 2n+2 {0) = X 2n+1 Vn G N. 

(3a) If a 0, then for any n > 



<9g 

lfn = 0, then ()26|) is va/zc? /or any a. 
(3b) If a > (respectively, a < 0), then 



° En (0) = 0. (26) 



d 2 E n/ . 1 / 4^(0) + 1 /■ \ \ 

r (°) = Tw^Tl 7777^ o., , - 8^(0) (27) 



<9g 2 8V6G(| - K(0)) V 

/or any even (respectively, odd ) n. 7/n = 0, iaen (}2Tj) zs valid for any a. 

(3c) //a = 0, t/ien (dE 2n+1 /dq) (0) < 0, (<9£ 2n+2 /<9g)(0) > 0, and (\(dE 2n+1 /dq)(0)\ = 
\(dE 2n+2 /dq)(0)\ Vn G N. 

(4) Ifq^O, then A n _i < £ n (g) < A n Vn G N. 

(5) lim £ n (g) = A n Vn G N. 

q— >oo 

Proof. Item (4) follows immediately from Proposition 3 and the definition of the functions 
E n . Formula (|19p shows that Q((; 0) = if and only if ( = A 2n+ i for some n G N; therefore 
items (2a)-(2d) follow from Theorem 2. Using the standard version of the implicit function 
theorem and the Proposition 3 again, we see that E n (q) are real analytic at q > 0. Moreover, 
item (3) of Theorem 2 implies that (i) E n (q) are real-analytic at q = for even n if a > 0, 
(ii) E n (q) are real-analytic at g = for odd n if a < 0, and (iii) -Eo(g) is real-analytic at 
g = for any a. In all these cases, the derivatives of E n can be found from the equations 

dQ dE n + <9Q _ Q 
d( dq dq 

dQ d 2 E n + cfQ fdK\ 2 + 2 & 2 Q_dK | d 2 Q _ Q 



<9£ <9g 2 <9£ 2 \ dq J d(dq dq dq 2 

Since (dQ/dq)(E; 0) = if E £ spec(H°), equation fl2BJ) follows from (|25|) in the consid- 
ered cases. In virtue of (J2HJ), the second derivative of E n is given by 

^OH-SfSVWo). (29, 



<9g 2 dq 2 \ d( 
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Substituting (J2H and <fifty in and using @ we get (|2Tj) . 

Now consider the singular case when -E n (0), n > 1, coincides with a point of the form 
\2m+\- In a neighborhood of the point (E n (0),0), introduce the function 

QaiCq) = — tt 

which is smooth with respect to (£, g) and analytic with respect to the first argument (. In 
a vicinity of (E n (0), 0) we have 

Q a (E n (q);q) = 0. (30) 

Further, 

«9Q a _ 1 <9Q r (1/2 - Q , , 

<9C r(i/2-c)SC J r2(i/2-C)' 1 J 

Since (dQ/d()((; 0) is a finite number at £ = A 2m +i and T(l/2 — £) has a pole at A 2m +i, 
the first term in (p!Tj) vanishes at the point (A2 m +i,0). The value of the function T'(l/2 — 
C)/r 2 (l/2 — at ( = A 2m +i is a nonzero finite number. Finally, Q(A 2m +i,0) = 0; thus 
dQa/dq vanishes at the point A2 m +i if and only if a = 0. Therefore, if a ^ 0, then each 
function E n (q) has an analytic continuation in a neighborhood of the point q = 0. Since 
q i — > Q a {^2m+u q) is an even function, we get easily (j2El)- 
Let now a = 0. Then 

d 2 Q = 1 d 2 Q 
d( 2 r(l/2 - C) <9C 2 

r(i/2-Q gg r-(i/2-c)r(i/2-o-2r 2 (i/2-c) 
r 2 (i/2-c)9C w r3(i/ 2 -o ' 1 1 

It is easy to see that the first and last terms in (f3*2*j) vanishes at the point (A2m+i> 0), whereas 
the second one does not. Therefore, d 2 Qo/d( 2 ^ at the point (A2 m +i,0), and E n (q) being 
solutions of ([30|h are continuous branches a two- valued analytic function in a vicinity of 
(A2 m +i,0). Obviously, at the point (A2 m +i, 0) the following relation is valid 

d 2 Q fdE n \ 2 n d 2 Q dE n dQ d 2 E n d 2 Q _ n 



<9£ 2 \ dq J d(dq dq d( dq 2 dq 2 
Since dQo/d( = at the considered point, we get the quadratic equation for dE n /dq: 

d 2 Q fdE n \ 2 d 2 Q 



(9C 2 \ dq J dq 

As a result, we complete the proof of items (1) and (3c). It remains to prove (5). Fix 
n G N and let e, < e < 1, is given. According to (J23J) we can choose q > such that 
<5(A n — e; q) — a < if q > q . Since Q(E n (q); q) — a = and the function i— > Q(E; q) 
increases in the interval \ n -\ < E < A n , we have E n (q) > \ n — e as q > q . Moreover, 
E n (q) < A n , and the proof is completed. I 

The structure of spec(if a (g)) given by Theorem 3 is presented in Table [T] The peculiar- 
ities of this table at q = can be understood from the point of view the symmetry group of 
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g = 




a > 


<Tl = {^2n(0) :ii6N} 

o"2 = {A 2 „+2 :i£N} 
0-3 = {A 2 n+i :i£N} 
a 4 = 


ax = {£ n (g) : iiGN} 
0-2 = {A n+ i : n G N} 

a 3 = 

a 4 = 


a = 


a x = {1/2} 

^2 = {A 2 n+2 :i£N} 

a 3 = 
o"4 = {A 2n +i :n£N} 


o x = {E n {q) : n e N} 
^2 = {A„+i :n£N} 

a 3 = 

a 4 = 


a < 


a t = {E 2n+1 (0) : neN}U{£ (0)} 
o"2 = {A 2n+2 : n e N} 
o"3 = {A 2n+ i :f»£N} 
<r 4 = 


ax = {£ n (g) : n G N} 
02 = {A n+ i : n G N} 

a 3 = 

a 4 = 



Table 1: The structure of spec(H a (q)) . 



the problem. It is well known that for a generic spherically symmetric potential V(r), the 
eigenvalues A of the operator H° = —A + V are parameterized by three quantum numbers: 
A = \ nr ,i,m, where n r (n r = 0, 1, . . .) is the so called principal (or total) quantum number; 
/ (I = 0, 1, . . .) is the orbital quantum number, and m (m — —I, — I + 1, ... ,1 — 1,1) is the 
magnetic quantum number. Each eigenvalue \ nr j^ m is degenerate with multiplicity 21 + 1, 
namely, \ nr j^ m = \ nr ,i,m> if rn,m' G {—/,—/ + 1, ... ,1 — 1,1}. This degeneracy is related 
to the invariance of H° with respect to the rotation group SO(3): eigensubspaces of H° 
carry an irreducible representation of this group. In general, X nr ,i,m ^ ^n' r ,i',m' if n r ^ n' r 
or / ^ V . The eigenvalues of an isotropic harmonic oscillator have an additional (so-called 
accidental) degeneracy: each eigensubspace L n is decomposed on the subspaces with 
angular momentum I = n, n — 2, . . . , (if n is even) or I = n,n — 2, . . . ,1 (if n is odd). 
This accidental degeneracy is related to the invariance of the Hamiltonian H° of an isotropic 
harmonic oscillator with respect to the group U(3). Indeed, 

3 3 
i=i 

where al- and a 3 - are standard creation and annihilation operators.^ Therefore, H° is invari- 
ant with respect to the transformation 

3 3 
3=1 3=1 

where (v,jk) is a unitary matrix. If q — 0, then H a (0) is a spherically symmetric perturbation 
of H° that violates the U(3)-symmetry. To prove this, we note that operators aja^ are gen- 
erators of the Lie group u(3). Therefore, if H a (0) is invariant with respect to the considered 
representation of U(3), we must have [H a (0), H°]—0. On the other hand it is easy to show 
that for C e C \ R the operator [R a ((), R°(()] has a nonzero integral kernel. 
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Since point perturbations can not change states with nonzero angular momentum / (see, 
e.g., IS]), the part cr 2 (at q = 0) may contain only even eigenvalues A 2n and we see this in 
Tabled Since all states from L n have the same parity (— l) n , the isotropic oscillator has no 
stationary states with a nonzero dipole momentum.^ On the other hand every eigensubspace 
of Hq(0) with eigenvalue from 04 have an eigenfunction with I = (this is the eigenfunction 
from item 4 of Theorem 1). Therefore point perturbations of an isotropic harmonic oscillator 
can lead to an appearance of eigenstates with nonzero dipole momentum. 

An alternative tool to understand the energy degeneracy of the three-dimensional isotro- 
pic oscillator gives the supersymmetry theoryP^"^"^We will not dwell here on this approach, 
nevertheless note that the analysis performed in the cited papers requires a modification in 
the s- channel only. 

The functions E n depend not only on the position parameter q, but also on the strength 
a; we will denote these dependencies as E n = E n (q,a). If E(q,a ) coincides with one of 
the numbers £ m , then in a vicinity of «o, the function a 1— »■ E n (q, a) is a continuous branch 
of the inverse function to E 1— > Q(E;q). It is already known from Proposition 1 that the 
following limits take place: 

lim E n (q;a) = X n , lim E n (q; a) = A n _i, 

a— >+oo a— *— 00 

where A_i = —00. Now we make more precise this behavior. From (|22j) we get the asymp- 
totics of the function Eo(q; a) for the fixed q > as a — > —00, 



E (q; a) = -167r V + \ + + O — , (33) 



q - 1 q (1 

4 + 128vr 2 a 2 + 

or in terms of the point perturbation of the initial operator (jlj 



%«) = -^ + ^ + ^ + 0(ffl, (34) 



where the scattering length l s tends to 0. Expression means that up to the infinitely 
small term 0{l 2 s ) the ground state of H a (q) equals to the ground state of the point per- 
turbation of the free Hamiltonian —h 2 A/2fi with the same scattering length l s shifted by 
the potential V(r) = fiQ 2 r 2 /2 at the point r = q. Equation (|34j) shows that at least for 
the isotropic harmonic oscillator its potential can be recovered from the dependence of the 
ground state of the point perturbation on the position of the potential support. It is reason- 
able to suppose that this is true for more general forms of the potential V; we consider this 
conjecture elsewhere. 

Now consider the behavior of E n (q; a) in a vicinity of the poles of Q((, q). We start with 
the general case q 7^ 0. Using (|2T|) we get as a — > ±00 



E,( q -,a) = Xi- ^ /2 ^l (2(n + l)^ n ( q 



+V2(g-' - q) H n (q/Vi) ff„+i (9/V2) + H?, +1 (9/V2)) a" 1 + 0( a - 2 ), (35) 
where A^ = A n and n > as a —>■ +00, and A~ = A n _i and n > 1 as a — > —00. 
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In the case q = 0, we are in position to give a compact form for more precise asymptotics 
of E n (q;a). Denote 



-*(2n)\\J \ 2^k(l + 2k)J 



(27r)S(2n)!! \(2tt) 
For eigenvalues with even indices we have 

\2n-1 for a < and n > 1, 



E 2 n(0-a) = { 



For the odd indices 



-16ttV + — ]—a- 2 + 0(a- 4 ) for a ->• -00 and n = 0, (36) 
12o7r 

A2n — A n (a) + 0(a~ 3 ) for a — > +00 and n > 0. 



A 2n +i for a > 0, 

^ 2 n+i(0;a) = { (37) 

— A n (a) + 0(a 3 ) for a — > —00. 

Formulas (|33j) - (|3T|) explain peculiarities in the plots of functions E n on Figures El and |U 
Note that in Eqs. (}3*3*|) - (}3T|) the remainder terms depend on n. 

The isotropic harmonic oscillator has an equidistant spectrum. After the perturbation 
by a zero-range potential, the distances between energy levels are changed and become 
dependent on the energy index n. This is important in the connection with the problem of 
the controlled modulation of the binding energy of the impurity center in quantum dots, that 
can be used to design nonlinear opto-electronic active elementsP The asymptotic formulas 
()34|) - (}3T|) give very accurate expressions for the excited energies in the most interesting 
case of a deep zero-range well (a — > —00) as well as for the case of a shallow well (a — > +00), 
which confirm numerical results from Ref. El Note also that Proposition 1 and Theorem 3 
imply a remarkable distinction between the excited energy for the ground state and that for 
the other ones: the energy E\(q; a) — E (q; a) can take an arbitrary value depending on q 
and a; on the other hand, energies A n — E n (q; a) and E n+ i(q; a) — A n (n > 1) are bounded by 
1. Since at fixed a, a —1, the function q \— > Ei(q; a) — E (q; a) is injective for moderate 
values of q, the position of an impurity in the quantum dot may be determined from the 
spectroscopy data. 

We show the plot of the energies E±(q; a) — E (q; a) and Ai — Ei(q; a) as functions of q 
and a on Fig. |S] and Fig. respectively. 

In conclusion we give the following remark. Let q : [0, 00) — > R be a smooth function 
obeying the conditions 

(HI) q > and the function r q(r) + r 2 /4 is nondecreasing, 
(H2) q'(r) < 0, 

and let kq and K\ are the first two eigenvalues of the operator H° + q = —A + r 2 /4 + q(r). 
It is proven in Ref. HHltliat kq/k\ < A0/A1, if q 7^ 0. Using Theorems 3 and A it is easy to 
construct smooth functions q with properties (HI) and 
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Figure 6: The exciting energy as a function of a for (a) q = 0, (b) q — 1, (c) q = 3. 

(H2a) q\r) > 0, 
such that Ko/ki > Ao/Ai. 
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